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ABSTRACT 



As microcircuit technology advances, there is an 
increased need for the dissipation of the heat which is 
generated. Extended surfaces are a useful tool in fulfilling 
the heat sink requirements for a microcircuit element. Heat 
transfer is very effective in the boiling mode and this thesis 
focuses on the analysis of a spine-shaped extended surface in 
a boiling liquid. Because the heat transfer in the spine and 
the propagation of signals on a transmission line are governed 
by identical differential equations, an analysis procedure 
based on the cascading of electrical transmission lines is 
used to predict the performance of the spine. The result of 
the analysis is a computer program that can assist the circuit 
designer in meeting any heat transfer requirements. 
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I . INTRODUCTION 

A major concern in the design of electrical equipment is 
the dissipation of heat inherent in the use of the equipment. 
Specific measures must be taken to transfer the heat from the 
circuit element (s) where it originates to the ultimate heat 
sink. This thesis investigates the performance of a spine that 
augments the transfer of heat from a highly dissipating 
circuit element to a boiling liquid. 

Extended surface technology has provided a means to 
transfer heat from microcircuit elements effectively and 
efficiently. There are many extended surface geometries that 
provide excellent heat transfer augmentation. However, this 
study will concentrate on the analysis of the spine geometry. 

The study focuses mainly on a computer program based on a 
cascade analysis deriving from the theory of transmission 
lines. The algorithm will enable the user to obtain a 
prediction of the surface temperature of an operating 
microcircuit element. This prediction is based on certain 
input parameters such as the spine geometry and thermal 
properties. This can potentially aid in the circuit design 
process where a heat sink must provide a certain overall heat 
transfer capability. If the rate of heat transfer is not 
adequate, other factors such as the material, surrounding 
temperature, heat flow mechanisms, and spine dimensions can be 
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adjusted and manipulated in order to arrive at desired 
dissipation rates. 

A primary application for a spine heat sink is in the 
thermal control of a microcircuit element. In an electric 
circuit, the power source generates the required energy so 
that the individual circuit elements may carry out their 
function. In some instances the energy required may be 
significant. A large amount of this electrical energy is 
converted to thermal energy. This thermal energy can create, 
because of the very small size of the microcircuit, a heat 
flux which approaches that of the surface of the sun 
(5xl0 6 W/m 2 ) . An extended surface such as a spine can be 
employed to augment the transfer of much of this heat so that 
the circuit will not overheat and can operate with the desired 
efficiency. The thrust of this thesis is to establish an 
analogy between the electrical transmission line and heat flow 
in a spine and apply the electrical theory to the analysis of 
a thermal system. It has already been proved that electrical 
transmission lines can be analyzed by subdividing a lossy line 
into segments and then cascading these segments using a 
transmission matrix. After the equations for the 
electro-thermal transmission line analogy are presented, it is 
theorized that the spine can be analyzed using a similar 
cascade procedure. A computer program is developed based on 
this cascade analysis concept for a "thermal transmission 
line", in this case the spine. 
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The electrical transmission line concepts are included 
solely to present enough information to allow for a comparison 
to the thermal transmission line. It is not intended to be a 
rederivation of the transmission line equations. 
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II. REVIEW OF EXTENDED SURFACES 



The earliest studies in the extended surface technology 
apparently occured in the year 1922. A paper by Harper and 
Brown (1922) reported on the mathematical analysis of heat 
transfer by conduction and convection occuring in a single 
extended surface. This extended surface became known simply as 
a fin. Some early experiments included coating metal rods with 
wax and observing the thermal properties as the bases of the 
rods were heated and the wax melted. The French mathematician 
Fourier (1822) also published mathematical analyses of the 
temperature variation of thin metal bars or rods, but these 
were not considered as fins or extended surface. Indeed, the 
work by Harper and Brown (1922) is considered to be the 
forerunner which has resulted in the advance of this 
significant subject area in the general field of heat 
transfer. From this modest beginning, the analysis and 
evaluation of extended surfaces has blossomed from a simple 
fin to arrays of surfaces where individual components are 
assembled into elaborate configurations. 

Some typical examples of extended surface are displayed in 
Figure 2.1. It can be observed that the more typical 
geometries are rectangular or spine shaped. Depending on the 
application, variations of these geometries can be formed to 
suit individual design needs. The extended surfaces shown are 
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Figure 2.1. Some typical examples of extended surfaces 

but a small sampling of those currently in use and/or under 
investigation. 

Specific terminology to be used is shown in Figure 2.2 for 
the spine of cylindrical profile which is sometimes referred 
to as the pin fin. Attention should be focused on what is 
meant by fin height. This is the distance between the fin base 
and the fin tip. The fin diameter is simply the diameter of 
the cylindrical spine and is assumed to be constant. It is 
noted that the origin of the height coordinate is at the fin 
tip with positive orientation in a direction toward the base. 
If the prime surface (also called the base surface) is at a 
higher temperature than the environment, the heat flow in the 
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Figure 2.2. A cylindrical spine 
[from Ref. 9] 

spine is from base to tip in a direction that is opposed to 
the positive sense of the coordinate system. 

It was determined by the early investigators in the 
extended surface technology that a set of assumptions were 
needed to make the mathematical analysis tractable. These have 
been known since 1945 as the Murray-Gardner (1945) assumptions 
and are taken from the papers by Murray (1938) and Gardner 
(1945) . The assumptions have considerable importance to 
ongoing analyses and research because their elimination, 
either one at a time or in any combination, provide a series 
of guidelines for subsequent investigators to follow. The 
following assumptions regarding extended surface analysis are 
taken directly from the Gardner (1945) paper. 

- The heat flow and temperature distribution throughout 
the fin are independent of time, i.e., the heat flow is 
steady. 

- The fin material is homogeneous and isotropic. 
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- There are no heat sources in the fin itself. 

- The heat flow to or from the fin surface at any point is 
directly proportional to the temperature difference 
between the surface at that point and the surrounding 
fluid. 

- The thermal conductivity of the fin is constant. 

- The heat transfer coefficient is the same over all the 
fin surface. 

- The temperature of the surrounding fluid is uniform. 

- The temperature of the base of the fin is uniform. 

- The thickness is so small compared to its height that 
temperature gradients normal to the surface may be 
neglected. 

- The heat transferred through the outermost edge of the 
fin is negligible compared to that passing through the 
sides . 

Kern (1950) added another assumption to Gardner's list: 

- The joint between the fin and the tube or prime surface 
is assumed to offer no bond resistance. 
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FIN EFFICIENCY 



The fin efficiency has been defined as the ratio of the 
heat transferred from the fin to the heat that would be 
transferred by the fin if its thermal conductivity were 
infinite (if the entire fin were to operate at the base 
temperature) . Fin efficiency is a function of the fin 
dimensions and thermal properties. Gardner (1945 1 ) defined a 
"fin effectiveness" as the ratio of the heat transferred into 
the base of the fin to the heat transferred through the same 
prime or base surface area if the fin were not present. He 
also determined a relationship to permit the conversion from 
fin efficiency to fin effectiveness. Although the concept of 
fin effectiveness has been in guestion, the theory is guite 
simple . 

Using 0 b = T b - T s as the base temperature excess above 
the surroundings in the heat transfer rate eguation, 

q = hS6 b W 

the surface, S, to be used is composed of fin surface S f and 
base or prime surface S b . Because of the temperature gradient 
within the fin, all of the fin surface does not operate at 0 b . 
The fin efficiency T) modifies the fin surface to allow the 
total surface area S to be represented by 
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S — Sf, + T]Sf 



( 2 ) 



so that for ease of computation all of the surface can be 
presumed to operate at the base temperature excess, 0 b . 
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HEAT TRANSFER IN THE BOILING MODE 



Boiling may occur when a surface is immersed in a liquid 
and is maintained at a temperature above the saturation 
temperature of the liquid. The heat flux will depend on the 
difference in temperature between the surface and the 
saturation temperature. Subcooled or local boiling occurs if 
the bulk temperature of the liquid is below the saturation 
temperature. If the liquid is sustained at the saturation 
temperature, the process is referred to as saturated, or bulk, 
boiling. 

The boiling curve is shown in Figure 4.1 where heat flux 
data from an electrically heated platinum wire submerged in 
water are plotted against temperature excess. The heat 
transfer data for boiling water was first obtained by Collier 
(1972) . The temperature excess in this case is the temperature 
difference between the heated wire and the surrounding liquid. 
In region I, the liquid near the heated surface is superheated 
slightly at the lower temperatures. The heated water 
subsequently rises to the surface where it evaporates. In 
region II bubbles begin to form on the surface of the wire and 
dissipate into the liquid as they break away from the surface. 
This region marks the beginning of the nucleate boiling 
regime. With a further increase in temperature excess, bubble 
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Temperature excess, 0, °C 



Figure 4.1 The boiling curve 

formation is more rapid. These bubbles break away and rise to 
the surface where they are dissipated. This is a 
characteristic of region III. 

Eventually, the bubble formation becomes so rapid that 
they totally envelop the heating surface. This prevents any 
inflow of liquid from taking their place. At this point the 
bubbles coalesce, forming a thin vapor film which blankets the 
surface. This vapor film acts as an insulator between the 
heated surface and the liquid. The heat must be conducted 
through this film to the liquid in order to effect the boiling 
process. The thermal resistance of this insulating film causes 
a dramatic reduction in heat flux. This phenomenon is known as 
transition to film boiling and is illustrated in region IV. 
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Note the negative slope of this region as compared to the 
positive slope of region III. Region IV is the transition 
region from nucleate boiling to stable film boiling and is 
unstable. 

Region V depicts stable film boiling. The surface 
temperatures required to maintain stable film boiling are 
high. Much of the heat lost in this region may be the result 
of thermal radiation. This region shows the return of the 
positive slope on the boiling curve. The curve in region V 
continues upward until the melting point of the heated surface 
is reached. 

The boiling process is unstable at point "a" on the 
boiling curve because a decrease in the heat flux results from 
just a small increase in temperature excess. Because the 
increased heat flux requirement is not met, operation further 
along the boiling curve results. Eventually, equilibrium can 
again be attained at point "b" in the film boiling region. It 
can be observed that an increase in temperature excess at this 
point yields a corresponding increase in heat flux. However, 
this is rarely seen in practice because this temperature 
usually exceeds the melting point of the material and burnout 
results. It is obvious that the preferred regions of the 
boiling curve for heat transfer are in regions II and III 
where nucleate boiling exists. 

Detailed photographs of the boiling modes using methanol 
as the liquid coolant are shown in Figure 4.2. They 
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Figure 4.2. 



Photographs of the 
[from Ref. 9] 



boiling modes 
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illustrate the energetic nucleate boiling regime. Transition 
boiling is seen in Figure 4.2b as the bubbles begin to 
coalesce at higher temperatures. Figure 4.2c shows film 
boiling and the thin vapor layer can be seen as well as the 
dramatic decrease in bubble formation. 

A. FIN OPTIMIZATION IN THE BOILING MODE 

The optimization of a cylindrical spine dissipating heat 
to a boiling liquid was thoroughly studied by Haley and 
Westwater (1966) . Their determinations of preferred heat flow, 
temperature gradient and spine dimensions were selected for a 
minimum spine volume rather than a minimum spine area profile. 
A constant thermal conductivity was assumed even though, 
realistically, a variation in exterior heat transfer 
coefficient and fin material thermal conductivity usually 
exists . 

It was found that several types of boiling existed 
simultaneously on the spine. The boiling forms most often 
encountered were nucleate and film boiling. This was the 
result of different heat transfer coefficients for various 
sections of the spine. A numerical solution for the heat 
transfer coefficient was derived using boiling curves, 
Simpson's rule and a fourth order Runge-Kutta approach. The 
optimum spine resulting from the previous determinations 
possesses a spadelike shape. 

The unusual spine geometry is justified when the 
distribution of the surface heat flux is considered. Heat 
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transfer is greater in the nucleate and transition boiling 
modes than it is in the natural convection and film boiling 
modes. To keep the surface area in the film boiling region to 
a minimum, the base of the spine has a small cross section. 
This enables the transition from film boiling to nucleate 
boiling to occur over a short distance. At the point where the 
heat transfer coefficient increases, the spine diameter is 
increased dramatically. This allows for a larger surface area 
to be exposed to the liquid operating in the transition and 
nucleate boiling regions. At the tip where the heat transfer 
coefficient decreases with lower temperatures, the spine 
tapers to a point. The spine with this shape makes optimum use 
of the most effective regions on the boiling curve. 

It was pointed out by Siman-Tov (1970) that the spade 
shaped spine had several undesirable attributes. Some major 
concerns were the difficulty and expense to manufacture such 
a design. Moreover, the spine would be difficult to weld and 
when attached to the prime surface, the small spine base could 
make the structure too fragile. Versions of this optimum spine 
were proposed that would include the basic design theory and 
also overcome some of these operational problems. Siman-Tov 
(1970) introduced a series of longitudinal fins protruding 
from a cylinder and a disk-type fin attached to a cylindrical 
spine in a "thumbtack" type of configuration. Cash, Klein, and 
Westwater (1971) concluded that two cones attached to a small 
cylindrical neck was a good approximation to the optimum 
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turnip-shaped spine. A computer was used to predict the heat 
transfer capabilities of this two-cone arrangement. As 
expected, this arrangement had greater heat transfer rates per 
unit volume than those calculated for a cylindrical spine. 
Experiments proved that the computer analysis was conservative 
and that the actual fins gave higher performance than the 
computed prediction. This, and the other works discussed, 
showed that the use of extended surfaces in the boiling mode 
can significantly increase the heat dissipation. 
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V. MATHEMATICS OF EXTENDED SURFACE ANALYSIS 



The heat transfer characteristics of individual fins are 
analyzed by a procedure summarized by Kern and Kraus (1972). 
The procedure begins with the differential equations relating 
the temperature excess and heat flow in each fin. The 
differential equations are solved using initial values for the 
temperature excess and heat flow. When many fins are 
physically coupled the problem then becomes one of solving a 
set of linear algebraic equations after the individual 
differential equations are integrated analytically. 

The procedure developed here treats each fin in the 
configuration as a lumped parameter matrix. Each fin or fin 
section can then be combined via matrix operations. Some 
improvements of this algorithm over existing techniques are as 
follows : 

- The differential equations for the individual fins are 
uncoupled. 

- The process which mathematically assembles the 
characteristics of the entire configuration from the 
characteristics of the individual fins is simple. It may 
even be done on a hand calculator and the effect of each 
parameter is easily isolated. 

- The performance of the individual fin can be measured by 
the ratio of the total heat dissipated to the 
temperature excess existing at the base of the fin. 

- While the analysis evaluates heat exchange operating 
under the Murray-Gardner considerations, some of these 
may be relaxed. The two Murray-Gardner assumptions which 
are relaxed in this procedure are as follows: 
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* There is no contact resistance between the base of 
the fin and the prime surface. 

* The heat transferred through the outermost edge of 
the fin (the fin tip) is negligible compared to 
that transferred through the lateral surfaces of 
the fin. 

A. THE HEAT FLOW IN A SPINE 

The heat transfer texts tell us that for the cylindrical 
spine displayed in Figure 2.2, the differential equation for 
the temperature excess, 0, above some datum level (the 
environmental temperature) will be 

dQ 2 2 

— -m 2 0 = 0 (3) 

where 

m = -Ja hTkd ( 4 ) 

Because Eq(3) is an ordinary second order differential 
equation, it has a general solution given by 

e{x) = C 1 e mx + C 2 e~ mx (5) 

According to Fourier's law, the heat flow will be 

q{x) = M^ = kAm^e™ - C 2 e~ mx ] ( 6 ) 

because the height coordinate orientation is opposed to the 
heat flow. 

The arbitrary constants C 1 and C 2 in Eqs(5) and (6) can be 
evaluated from the initial conditions 
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6{x -0)- d a -C 1 +C 2 ( 7 ) 

and 

q(x = 0) = q a = hAm[C 1 -C 2 ] ( 8 ) 

When these are inserted into Eqs(5) and (6), two expressions 
containing C,, and C 2 are obtained 



Ci +C 2 - 6 a 



and 



Ci-C 2 = -^_ 

hAm 

The addition of Eqs(9) and (10) gives 



2Ci — 6 n + 






a hAm 



= e a+ Y 0 <la 



where Y - 1/hAm = (7r 2 hkd 3 ) h /2 



(9) 



( 10 ) 



( 11 ) 



Solving for C, and C 2 yields 



and 



Cl =|[e.+Vrf.] 



( 12 ) 



C 2 = e a - Y 0 q a 



(13) 



When these constants are substituted into Eqs(5) and (6), the 
expressions for the temperature excess and heat flow can be 
written as hyberbolic sines and cosines. The temperature 
excess and heat flow at any point "x" will be 
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0(x) = O a cosh mx + Y 0 (j a sinh mx 



(14) 



and 



( ](x ) = hAmd g sinh mx + hAmY 0 cosh mx 



(15) 



In particular, when x=b a linear transformation, mapping 
conditions from the spine tip to the spine base, can be 
proposed and the matrix of this linear transformation is 



B. THE TRANSMISSION LINE ANALOGY 

With the use of transmission line concepts, the 
distributed heat flow pattern of a spine can be represented by 
an analogous distributed electrical network. To analyze this 
system, the transmission line can be represented as a series 
of lumped parameter "tee" sections. Consider a line of length, 
L, with a length coordinate, x, having its origin at the 
receiving end and its positive orientation toward the sending 
end. Such a line is shown in Figure 5.1a where the subscripts, 
S and R, refer to the sending and the receiving ends 
respectively. An isolated differential section of transmission 
line is displayed in Figure 5.1b. 



T = 



cosh/nb Z 0 sinh/nb 

Y 0 sinh/nb coshmb 



(16) 
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« 7 ) 




(b) 

Figure 5.1. The electrical transmission line 

[ f r om Ref. 9 ] 



The two Kirchoff laws may be used to obtain the so-called 
'telegraph equations' which relate the voltage and the current 
at any time and at any point along the line. These telegraph 
equations are 

^ = LC ^ + ( RC + LG )^ +RGv (17) 

and 



dx 2 LC li 2 +( - RC+LG ^Tt +Ra 



(18) 
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Consider a lossy line in which L=0 and C=0. In this event 



Eqs ( 17 ) and (18) become 



and 



d 2 v 

dx 2 



-ccv = 



0 




0 



(19) 



( 20 ) 



where the partial derivative is no longer needed because 
neither v nor i depend on time. These ordinary second order 
differential equations have general solutions 



v(x) = V 1 e ax + V 2 e~ ax ( 21 ) 

and 

i(x) = l 1 e ax + l 2 e~ ax ( 22 ) 

Use of the initial conditions at the receiving end of the 
line where x=0 gives 



®(°) = Vr 

and 

l (°) = 



At any point x along the line, 



r= Ri 

dx 



and 




(23) 



(24) 



(25) 



( 26 ) 
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and these, coupled with Eqs(23) and (24), yield the total 
solution for the voltage and current at any point along the 
line 



i>(x) = V R cosh ax + Z 0 1 R sinh ax 



(27) 



and 



i(x ) = YqVr sinh ax + I R cosh ax 



(28) 



where Z is the characteristic impedance of the line 



Z n =VRG 



(29) 



and where Y 0 is the characteristic admittance of the line 



Yn = TT-= 1 



0 z n Vrg 



(30) 



The form of Eqs(27) and (28) suggests a linear 
transformation from conditions of voltage and current at the 
receiving end to any point along the line. 



'v(x)' 




cosh ax Z 0 sinh ax' 


~Vr 


i(x)_ 




Y q sinh ax cosh ax 


Jr. 



(31) 

In the case at the sending end where x=b this transformation 
can be written as 



( 32 ) 



'Vs' 




cosh ab 


Z 0 sinh ab * 


Vr 


Js. 




Y 0 sinh ab 


coshafe 


Jr. 
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A comparison of Eqs(16) and (32) shows the similarities. 
The linear transformation matrix for heat flow and temperature 
excess is almost identical to that for voltage and current in 
a transmission line. The only difference is that the value « 
in the electrical network is replaced by m in the thermal 
network. These parameters contain constants obtained from 
solutions to the respective differential equations. Therefore, 
there is an obvious analogy between the thermal transmission 
line and the electrical transmission line. 

Utilizing electrical two-port theory, the entire 
transmission line can be segmented into a series of minute yet 
finite sections. These sections can then be analyzed 
individually with the results of each assembled via a cascade 
procedure. This results in the treatment of each section as a 
two-port system indicated by Figure 5.2. The output of one 
section is the input to the next. This method has proven to 
be useful and accurate in the analysis of electrical 
transmission lines. 



q b w 
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^ — r~ 
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Figure 5.2. A two-port representation 
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It is logical that a similar cascade procedure may be 
successful in analyzing a fin or spine which is now referred 
to as a thermal transmission line. The linear transformation 
matrices are equivalent when « is replaced by m. In the 
thermal case, heat flow q and temperature excess 0 are the 
parameters which are "transformed" in each individual segment. 
Certain extended surface geometries, specifically the spine, 
allow for direct application of the cascade procedure for 
analysis. Two cascaded subfins are shown in Figure 5.3. A 
computer program was developed to test this procedure for a 
spine. 




Figure 5.3. Two cascaded subfins 
[from Ref. 9] 
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C. THE ALGORITHM FOR THE CASCADE PROCEDURE 

As mentioned earlier, a cascade algorithm was developed to 
analyze the heat transfer properties of a spine. The spine is 
divided into segments, or subfins, and each section is treated 
as an independent two-port system. The computations begin at 
the tip and conclude at the section nearest to the base 
surface. The inputs to the program are: 

- Temperature at the fin tip 

- Temperature of the surrounding liquid (boiling liquid) 

- Thermal conductivity of the spine material 

- Diameter of the spine 

- Height of the spine 

A flowchart for the cascade algorithm is shown in Figure 5.4. 
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Figure 5.4. Flowchart for the cascade algorithm 
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The following are the significant steps in the cascade 
algorithm: 

1) Assume the tip temperature, T a 

This value is entered by the user and must be 
greater than the saturation temperature. 

2) Compute the spine tip temperature excess, 0 a 

The temperature excess is found by subtracting 
the surrounding temperature from the tip 
temperature, 

0. - T. - V 

3) Find the heat-transfer coeficient, h 

The heat-transfer coefficient is dependent on the 
temperature excess, 0 a determined in step 2. The 
following are based on data given by 
Collier (1972) . 

• If 0 a is between^ 0.0 and 7.0 degrees, 

h = 1300. O(0 a )° , 

and the subfin is in the natural convection mode. 

• If 0 a is between 7.0 and 25.0 degrees, 
h = 76 . 8 (0 a ) 2 ' , 

and the subfin is in the pool-nucleate boiling 
mode. 

• If 0 a is between 2,5.0 and 280.0 degrees, 
h - 7.5 X 107(6.), 

and the subfin is in the transition boiling mode. 

• If 0 a is greater than 280.0 degrees, 

h = 8OO.O(0 a ) , 

and the subfin is in the film boiling mode. 



4) Compute the heat flow, q a 



The heat flow. q a , is derived from the relation, 
q a = h(7r/4)d 2 0 a . 
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5) Assume T b and calculate 0 b 

T b is given an initial value of T a 
0 b is then equal to 0 a 

6) Compute the average temperature for the subfin, T AV 

This is found by averaging T a and T b 
t av = ( T a +T b )/ 2 

7) Compute the average tmperature excess, 0 AV 

The average temperature excess is found using T AV 
from step 6, 0 AV = T AV - T s 

8) Compute the average heat-transfer coefficient, h AV 

The average heat-transfer coefficient is a 
function of 0 AV from step 7. As shown in step 3, 
the equation for calculating h AV depends on which 
boiling mode the subfin is experiencing. To find 
h AV , the equations in step 3 are used with 0 AV 
replacing 0 a . 

9) Calculate mb and Y 0 

This value is found using the relation 

mb = V 4h AV /kd (L/N) where L is the length of the 

entire spine and N is the number of subfins. 

The subfin admittance, Y 0 is found using the 

relation Y 0 = / 7r 2 hkd 2 /2 

10) Find the values for the transmission matrix, ' A,B,C,D' 



V 




A 


Ii~ 




A 




C 


D 


, 

a? 

1 



where 

A = cosh mb 
B = 1/Y 0 sinh mb 
C = Y 0 sinh mb 
D = cosh mb 

11) Update 0 b = 0^^ and compute q b 

The temperature excess, 0,^ is updated and heat 
flow, q b is calculated using the transformation 
matrix from step 10. 
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12) Compare the temperature excesses, 0 b and 0^^ 

The updated temperature excess, 0^^, from step 11 
is compared to the temperature excess, 0 b , from 
step 5. If the values are equal or within a 
margin of error, 0^^ becomes 0 a and q b becomes q a 
for the next subfin. The algorithm then returns 
to step 6 and continues. If the values are not 
equal, 0 b is given the value 0^^. A new value for 
T b is then computed and the algorithm returns to 
step 6 and continues. The last calculation in the 
subfin under analysis is qb/0 b , the input 
admittance, Y in . 



Y in = V 0 b = [C + D(q a /0 a )]/[A + B(q a /0 a ) ] 
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VI. THE CASCADE PROGRAM 



As previously indicated, the cascade program is based on 
a thermal transmission line model. The program that was 
developed is a menu-driven, user-friendly software package 
which will assist the user in analyzing heat transfer 
requirements. It is intended that the program be used as a 
tool to ascertain whether or not a spine is a suitable heat 
sink to meet the circuit design needs. 

The results of the program address two general 
design considerations. One is the adequacy of the heat 
transfer and the other is the compatability of the spine with 
regard to size limitations. Given the miniaturization of 
modern microcircuit packages, it is obvious that an effective 
heat sink may have strict size limitations. This chapter 
considers these aspects of the cascade program. Several sample 
problems are provided in order to give an assessment of the 
program's attributes. 

A. FEATURES 

The program is menu driven providing the user with a 
series of parameter inputs. Upon entering the program, the 
user will be provided with the option of viewing an optional 
overview. First-time users should consider reading this 
overview. The program is designed to operate with either upper 
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or lower-case letters. All units used in the algorithm are SI 
units . 

Following the overview option, the user is asked to enter 
the geometric parameters for the spine under analysis. The 
first input is the spine height. It is recommended that a 
height of approximately 25 mm be entered initially for base 
temperatures between 115°C and 200°C. The algorithm adjusts the 
spine height to the minimum required in order to meet the 
desired base temperature. The user then inputs the spine 
diameter which is not adjusted by the program. 

Spines may be manufactured from an assortment of 
materials, each with its own thermal conductivity. The program 
provides the option to test spines of different materials by 
entering the appropriate values. 

The user is then prompted to input the saturation 
temperature of the fluid. This value should be at least 100°C 
because the spine is analyzed in the boiling mode. After 
selecting the saturation temperature, the user must specify 
the desired base temperature. This is the surface temperature 
of the microcircuit element. The final entry is the desired 
spine tip temperature. The program may adjust this value 
slightly in order to meet the specified base temperature. 

With all of the data entered by the user, the program 
calculates key parameters concerning the spine's operational 
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characteristics. The spine base temperature is predicted 
within 0 . 1°C of the user's desired base temperature. The tip 
temperature may be adjusted slightly by the program and the 
temperature which is predicted is presented as output for 
observation by the user. The heat flow rate in watts is also 
provided as output. 

The spine is operating in the boiling mode and the user 
may wish to ascertain the boiling regimes that are predicted. 
The algorithm has the capability to distinguish the percent of 
the spine operating in natural convection, nucleate boiling, 
and transition boiling. The film boiling mode is not 
considered in the program because microcircuit applications do 
not require film boiling as a means to dissipate heat. 
However, a simple modification to the program will permit film 
boiling to be analyzed. 

B . SAMPLE PROBLEMS 

The following terminal session is a typical example of the 
cascade program and its capabilities. Six sets of parameters 
are tested and the respective outputs can be observed. 

»»»> RUN #1: 

Enter the height of the spine in millimeters: 20 

Enter the diameter of the spine in millimeters: 2 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 
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Enter the temperature of base surface in deg C: 110 
Enter the desired tip temperature in deg C: 102 



You have entered the following parameters: 



Spine height (mm) : 20.0 
Spine diameter (mm) : 2.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 110.0 
Tip temperature: 102.0 



Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



The computed spine base temperature is (deg C) : 109.8 



The spine tip temperature is (deg C) : 101.8 
The heat dissipated by the spine is (watts): 1.2 
The revised spine height is (mm): 19.0 
The spine diameter (user input) is (mm): 2.0 
Percentage of spine in natural convection is: 89.0 
Percentage of spine in nucleate boiling is: 11.0 
Percentage of spine in transition boiling is: 0.0 



Do you want to enter another data set? 

Enter "Y" for yes, "N" for no. 

Y 

»»>» RUN #2: 

Enter the height of the spine in millimeters: 20 
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Enter the diameter of the spine in millimeters: 2 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 

Enter the temperature of base surface in deg C: 115 
Enter the desired tip temperature in deg C: 102 

You have entered the following parameters: 



Spine height (mm) : 20.0 
Spine diameter (mm) : 2.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 115.0 
Tip temperature: 102.0 



Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



The computed spine base temperature is (deg C) : 115.0 
The spine tip temperature is (deg C) : 101.9 
The heat dissipated by the spine is (watts): 2.8 
The revised spine height is (mm): 20.0 
The spine diameter (user input) is (mm): 2.0 
Percentage of spine in natural convection is: 81.0 
Percentage of spine in nucleate boiling is: 19.0 
Percentage of spine in transition boiling is: 0.0 



Do you want to enter another data set? 
Enter "Y" for yes, "N" for no. 

Y 
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RUN #3: 



Enter the height of the spine in millimeters: 20 

Enter the diameter of the spine in millimeters: 2 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 

Enter the temperature of base surface in deg C: 125 
Enter the desired tip temperature in deg C: 102 



You have entered the following parameters: 



Spine height(mm): 20.0 
Spine diameter (mm) : 2.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 125.0 
Tip temperature: 102.0 



Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



The computed spine base temperature is (deg C) : 125.1 
The spine tip temperature is (deg C) : 102.0 
The heat dissipated by the spine is (watts): 7.8 
The revised spine height is (mm): 21.0 
The spine diameter (user input) is (mm): 2.0 
Percentage of spine in natural convection is: 75.0 
Percentage of spine in nucleate boiling is: 25.0 
Percentage of spine in transition boiling is: 0.0 
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Do you want to enter another data set? 
Enter "Y" for yes, "N" for no. 

Y 



>»»» RUN #4: 

Enter the height of the spine in millimeters: 30 

Enter the diameter of the spine in millimeters: 2 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 

Enter the temperature of base surface in deg C: 125 
Enter the desired tip temperature in deg C: 102 



You have entered the following parameters: 



Spine height (mm) : 30.0 
Spine diameter (mm) : 2.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 125.0 
Tip temperature: 102.0 



Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



The computed spine base temperature is (deg C) : 125.1 
The spine tip temperature is (deg C) : 102.0 
The heat dissipated by the spine is (watts): 7.8 
The revised spine height is (mm): 21.0 
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The spine diameter (user input) is (mm): 2. 
Percentage of spine in natural convection is: 75. 
Percentage of spine in nucleate boiling is: 25. 
Percentage of spine in transition boiling is: 0. 



Do you want to enter another data set? 

Enter "Y" for yes, "N" for no. 

Y 

»»»> RUN #5: 

Enter the height of the spine in millimeters: 20 

Enter the diameter of the spine in millimeters: 3 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 

Enter the temperature of base surface in deg C: 115 
Enter the desired tip temperature in deg C: 102 

You have entered the following parameters: 



Spine height(mm): 20.0 
Spine diameter (mm) : 3.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 115.0 
Tip temperature: 102.0 



Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



0 

0 

0 

0 
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The computed spine base temperature is (deg C) : 114.8 



The spine tip temperature is (deg C) : 102. 
The heat dissipated by the spine is (watts) : 5. 
The revised spine height is (mm) : 25. 
The spine diameter (user input) is (mm): 3. 
Percentage of spine in natural convection is: 81. 
Percentage of spine in nucleate boiling is: 19. 
Percentage of spine in transition boiling is: 0. 



Do you want to enter another data set? 

Enter "Y" for yes, "N" for no. 

Y 

»»»> RUN #6: 

Enter the height of the spine in millimeters: 20 

Enter the diameter of the spine in millimeters: 3 

Enter the thermal conductivity of the spine material 
in watts per meter per deg C: 202 

Enter the saturation temperature of the 
fluid in deg C: 100 

Enter the temperature of base surface in deg C: 125 
Enter the desired tip temperature in deg C: 102 

You have entered the following parameters: 



Spine height (mm) : 20.0 
Spine diameter (mm) : 3.0 
Thermal conductivity: 202.0 
Saturation temperature: 100.0 
Base temperature: 125.0 
Tip temperature: 102.0 



0 

0 

0 

0 

0 

0 

0 
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Do you want to make any changes? : 
Enter "Y" for yes, "N" for no. 

N 



The computed spine base temperature is (deg C) : 124.8 
The spine tip temperature is (deg C) : 102.0 
The heat dissipated by the spine is (watts): 14.0 
The revised spine height is (mm): 27.0 
The spine diameter (user input) is (mm): 3.0 
Percentage of spine in natural convection is: 76.0 
Percentage of spine in nucleate boiling is: 24.0 
Percentage of spine in transition boiling is: .0 



Do you want to enter another data set? 

Enter "Y" for yes, "N" for no. 

N 

Stop - Program terminated. 

It is observed that in all of the sample tests, the 
desired base temperature was attained with an accuracy of 
0.2°C. The spine height was adjusted by the computer code as 
necessary to obtain the ideal geometry for the desired heat 
dissipation. In some cases the spine tip temperature was 
adjusted slightly to assist in reaching the desired base 
temperature. 

Also apparent is the increase in heat flux in the spine as 
the base temperature increased. As expected, the heat flow 
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rate in the 3 mm diameter spine is greater than the heat flow 
in the 2 mm diameter spine for the same base temperature. 

The program output also shows the boiling modes which are 
predicted for the spine being tested. It is observed that, for 
the sample problems, most of the spine was in natural 
convection. A manipulation of the input parameters will 
permit a larger portion of the spine to operate in the 
nucleate boiling mode. 

The sample runs give an example of how the cascade program 
can assist a circuit designer in analyzing specific heat 
transfer requirements. The simple entries and quick results 
allow the user to test many spine geometries for an assortment 
of heat transfer needs. 
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VII 



CONCLUSION 



The purpose of this thesis is to develop a cascade 
algorithm to analyze a spine in the boiling mode. This cascade 
approach is based on an analogy between the heat flow 
equations in the spine and the equations of the electrical 
transmission line. 

The current version of the algorithm shows that cascading 
a series of subfins is an accurate and efficient method for 
analyzing an extended surface such as a spine. The results 
show that analyses can be made for a spine heat sink in a 
boiling environment without the need for time-consuming manual 
calculations . 

The expansion possibilities for the cascade program are 
significant. Presently, the algorithm is modeled for water as 
the boiling liquid. With some modifications, a user could be 
able to input data for other boiling liquids and the spine 
would be analyzed accordingly. Other potential features 
include the addition of graphics capabilities and the 
introduction of user-friendly peripherals such as mouse 
drivers. 

The possible applications for the cascade program are 
substantial. As the size of microcircuits decreases and the 
power requirements increase, the need for efficient and 
compact heat sinks is great. There is a real need for a tool 
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that allows the designer to accurately perform thermal 
analysis of microcircuits during the design process. The 
cascade program to analyze a spine is the first step in 
enabling the designer to effectively determine a heat sink 
that meets specified requirements. 
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APPENDIX A 
NOMENCLATURE 

, 2 
A Cross sectional area, (m ) 

A An element of the transmission parameter matrix, 

dimensionless 

a Coordinate of fin tip, (m) 

B An element of the transmission parameter matrix, (°C/W) 

b Coordinate of fin base, (m) 

C An element of the transmission parameter matrix, (W/°C) 

D An element of the transmission parameter matrix, 

dimensionless 

d Diameter, (m) 

h Heat transfer coefficient, (W/m 2 - °C) 

k Thermal conductivity, (W/m - °C) 

L Fin length, (m) 

m Fin performance factor, (m 1 ) 

q Heat flow, (W) 

2 

S Surface area, (m ) 

T Temperature, (°C) 

Y Admittance, (W/°C) 

Z Impedance, (°C/W) 

Matrix 

[T] Transmission matrix 
Greek 

. . . . -1 

« Transmission line attenuation, (m ) 

r) Fin efficiency, dimensionless 

0 Temperature difference or temperature excess, (°C) 
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Subscripts 
a Designates 

b Designates 

f Designates 

s Designates 



fin tip 
fin base 
fin 

surroundings 
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APPENDIX B 



PROGRAM LISTING 



The following is a listing of the cascade program used to 
analyze a spine operating in the boiling mode. It consists of 
main program followed by three subroutines. The key variables 
used in the program are given as well as comments throughout 



the algorithm. 



C 

C 

C 

C 

C 

C 
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C 
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C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c*** 



TITLE: 
AUTHOR : 
DATE: 
COMPILER: 
LINKER: 



CASCADED SPINE ANALYSIS 
LT JOHN W LARKIN 
03 MAR 1992 
MICROSOFT VERSION 3.2 
MICROSOFT VERSION 3.2 



MAIN PROGRAM 



PURPOSE: 

THIS PROGRAM PREDICTS THE PERFORMANCE OF A SPINE IN 
THE BOILING MODE 



DESCRIPTION OF MAIN PARAMETERS: 

ELINP — SPINE HEIGHT (mm) 

EL — SPINE HEIGHT (m) 

DIAM — SPINE DIAMETER (mm) 

DEE — SPINE DIAMETER (m) 

CONDTH — CONDUCTIVITY OF THE SPINE MATERIAL 
TSAT — SATURATION TEMPERATURE OF THE LIQUID 
TBASE — DESIRED BASE TEMPERATURE 
TBCALC — PREDICTED BASE TEMPERATURE 
TTIP — DESIRED TIP TEMPERATURE 
THA — TIP TEMPERATURE EXCESS 

THAV — AVERAGE TEMPERATURE EXCESS FOR SUBFIN 
THBN , THBO — BASE TEMPERATURE EXCESS 
TIPSUR — AREA OF SPINE CROSS SECTION (m A 2) 

H — HEAT TRANSFER COEFICIENT 
QA — HEAT FLOW AT THE TIP 
QB — HEAT FLOW AT THE BASE 
A — 'A' TRANSMISSION PARAMETER 
B — 'B' TRANSMISSION PARAMETER 
C — / C / TRANSMISSION PARAMETER 
D — 'D' TRANSMISSION PARAMETER 
EMB — FIN PERFORMANCE FACTOR (m) X FIN BASE 
COORDINATE (b) 

YO — ADMITTANCE 
ZO — IMPEDANCE 

PERNAT — PERCENTAGE OF SPINE IN NATURAL CONVECTION 

PERNUC — PERCENTAGE OF SPINE IN NUCLEATE BOILING 

PERTRN — PERCENTAGE OF SPINE IN TRANSITION BOILING 



*********** ***BEGINNING 
DIMENSION YIN (100) 



OF MAIN PROGRAM***************** 
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non ooooo oooo ooo 



CHARACTER* 3 CHNG 



4000 CALL CLS 

DISPLAY A PROGRAM OVERVIEW IF DESIRED 

4001 WRITE (IOT, 3008) 

READ ( IN , 2 8 4 0 ) ANS 



IF (ANS 


.EQ. ' Y' 


.OR. 


ANS 


.EQ. 


'y') THEN 


CALL OPTION 










GO TO 


4002 










ENDIF 
IF (ANS 


.EQ. 'N' 


.OR. 


ANS 


• EQ. 


'n') GO TO 4002 



GO TO 4001 



READ THE INITIAL PARAMETERS AND CALCULATE THE FIN 
TIP CONDITIONS 

4002 CALL CLS 

WRITE (IOT, 2803) 

READ (IN, 2802) ELINP 
WRITE (IOT, 2804) 

READ ( IN , 2802 ) DIAM 
WRITE (IOT, 2805) 

WRITE (IOT, 2811) 

READ (IN, 2802) CONDTH 
WRITE (IN, 2806) 

READ(IN, 2802) TSAT 
WRITE (IOT, 2807) 

READ (IN, 2802) TBASE 
WRITE (IOT, 3006) 

READ (IN, 2802) TTIP 



CORRECT THE INPUT VALUES IF THEY WERE ENTERED 
INCORRECTLY 



501 CALL CLS 

WRITE (*,3000) ELINP , DIAM , CONDTH , TSAT , TBASE , TTIP 
READ (*,3003) CHNG 

IF (CHNG .EQ. 'Y'.OR. CHNG .EQ. 'y') THEN 

502 CALL CORRECT 

RE AD (*,3004) INPUT 

IF (INPUT .EQ. 1 .OR. INPUT .EQ. 2 .OR. INPUT .EQ. 3 
+ .OR. INPUT .EQ. 4 .OR. INPUT .EQ. 5 .OR. INPUT 

+ .EQ. 6) THEN 

CALL UPDATE ( INPUT , ELINP , DIAM , CONDTH , TSAT , TBASE , TTIP) 
GO TO 501 
ELSE 

GO TO 502 
ENDIF 
ENDIF 

IF (CHNG .EQ. 'N' .OR. CHNG .EQ. 'n') GO TO 500 
GO TO 501 



SET THE INITIAL VALUES 



500 EL = ELINP/1000 . 
DEE = DIAM/ 1000 . 
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NAT = 0. 

NUC = 0. 

NTRANS = 0. 
PERNAT = 0. 
PERNUC = 0. 
PERTRN = 0. 
NCOUNT = 1 
FLAG = 0 
FLAGA = 0 
FLAGB = 0 
CHKR = 0. 

VALUE = 0. 

ADD = 0 

PI = 3.14159265 
BEE = EL/100. 



FIRST, OBTAIN THE TIP CONDITIONS 



THA = TTIP - TSAT 
10 THAI = THA 

TIPSUR = PI* (DEE**2. )/4. 

IF (THA. GT. 7. ) GO TO 60 
H = 1300* (THA** . 25) * ( (1 ./DIAM) ** . 25) 
MODE = 1 
GO TO 80 

60 IF (THA. GT . 40 . ) GO TO 70 
H = 76.8* (THA**2. ) 

MODE = 2 
GO TO 80 

70 H = 7.5*(10.**8.)/( THA* * 3 . ) 

MODE = 3 

80 QA = H*TIPSUR*THA 
N = 1 

THE TIP PARAMETERS HAVE BEEN CALCULATED 
PROCEED TO THE 100 SUBFINS 

THREF = THA 
90 THBO = THA 
95 THAV = (THBO + THA)/2. 

I F ( THAV . GT . 7 . ) GO TO 120 
H = 1300* (THAV**. 25) *( (1 ./DIAM) **. 25) 
MODE = 1 
GO TO 140 

120 I F ( THAV . GT .40.) GO TO 130 
H = 7 6 . 8 * ( THAV * * 2 . ) 

MODE = 2 
GO TO 140 

130 H = 7.5*(10.**8.)/ (THAV* *3 . ) 

MODE = 3 

140 EMB = SQRT ( 4 . *H/ (CONDTH*DEE) ) *BEE 

YO = SQRT (PI*PI*H*CONDTH*DEE*DEE*DEE) /2 
ZO = l./YO 

EPLUS = EXP (EMB) 

EMINUS = 1/ EPLUS 
A = (EPLUS + EMINUS )/2. 

D = A 

DUMMY = (EPLUS - EMINUS )/2. 

B = ZO*Dl IMY 
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no oooo non 



C = YO* DUMMY 

WE NOW HAVE THE T MATRIX FOR EACH SUBFIN 

THBN = A*THA + B*QA 
ZETZ = ABS (THBN - THREF) 

IF ( ZETZ .LE. .01) GO TO 175 
ADD = ADD + 1 
THBO = THBN 
N = N + 1 

IF (N . GT. 50 . AND. ZETZ . LE .6.5) GO TO 175 
IF (N.GT.50 .AND. ZETZ .GT. 6.5) THEN 
EL = EL - .01 
BEE = EL/100 
THA = THAI 
GO TO 10 
ENDIF 

THREF = THBN 
GO TO 95 

175 QB = C*THA 4- D*QA 
YIN (N) = QB/THBN 
IF (NCOUNT .EQ. 100) GO TO 300 
NCOUNT = NCOUNT + 1 
N = 1 
QA = QB 
THA = THBN 



IF 


(MODE . GT . 1) 


GO 


TO 


182 


NAT 


= NAT + 1. 








GO 


TO 90 








182 IF 


(MODE . GT. 2 ) 


GO 


TO 


184 



NUC = NUC + 1. 

GO TO 90 

184 NTRANS = NTRANS + 1. 
GO TO 90 



COMPARE THE PREDICTED BASE TEMPERATURE TO THE DESIRED 
BASE TEMPERATURE - ADJUST THE SPINE HEIGHT AS NECESARY 



300 TBCALC = THBN + TSAT 
ERR = TBCALC - TBASE 
ERRABS = ABS (ERR) 

IF (FLAGA .NE. 0 .OR. FLAGB .NE.O ) GO TO 305 
IF (ERRABS .LE. 3) GO TO 305 
IF (ERR .GT. 3 .AND. FLAG .EQ. 0) THEN 
EL = EL - .001 
BEE = EL/100. 

ENDIF 

IF (ERR .GT. 3 .AND. FLAG .NE. 0) THEN 
VALUE = VALUE + 1. 

EL = EL - .001/VALUE 
BEE = EL/100 
ENDIF 

IF (ERR .LT. 3) THEN 
EL = EL + .001 
BEE = EL/100. 

FLAG = 1 
ENDIF 

THA = THAI 

RESET VALUES FOR THE NEXT LOOP 

A Q 



non ooo ooooo 



c 



306 NCOUNT = 0 
ADD = 0 
NAT = 0. 

NUC = 0. 
NTRANS = 0. 
GO TO 10 



IF THE COMPUTED BASE TEMPERATURE IS WITHIN 
OF THE DESIRED BASE TEMPERATURE, ADJUST THE 
TEMPERATURE UNTIL THE ERROR IS LESS THAN 0. 



305 IF (ERRABS .LE. .186) GO TO 310 
IF (ERR .GT. .186) THEN 
IF ( FLAGB .NE. 0) THEN 
CHKR = CHKR+.15 
ENDIF 

THA = THAI - .05/ CHKR 
THAI = THA 
FLAGA = FLAG A + 1 
IF (FLAGA .GT. 100) GO TO 555 

GO TO 306 
ENDIF 

IF (ERR .LT. .186) THEN 
IF (FLAGA. GT.l) THEN 
CHKR = CHKR+.15 
ENDIF 

THA = THAI + . 05/ CHKR 
THAI = THA 
FLAGB = 1 
GO TO 306 
ENDIF 

DETERMINE THE SPINE'S BOILING REGIMES 



310 PERNAT = (NAT/100. ) *100. 
PERNUC = (NUC/100. ) *100. 
PERTRN = (NTRANS/100. ) *100. 



OUTPUT THE RESULTS 



CALL CLS 
WRITE (IOT, 2808) 
WRITE (IOT, 2809) 
WRITE (IOT, 3007) 
WRITE (IOT, 2809) 
WRITE (IOT, 2810) 
WRITE (IOT, 2809) 
WRITE (IOT, 2813) 
WRITE (IOT, 2809) 
WRITE (IOT, 2819) 
WRITE (IOT, 2809) 
WRITE (IOT, 2815) 
WRITE (IOT, 2816) 
WRITE (IOT, 2817) 
WRITE (IOT, 2816) 
WRITE (IOT, 2818) 
WRITE (IOT, 2816) 



TBCALC 
THAI + TSAT 
QB 

EL*1000 . 

DIAM 

PERNAT 

PERNUC 

PERTRN 



DEGREES C 
TIP 

DEGREES C 
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520 WRITE (IOT, 3005) 

READ(IN, 2840) ANS 

IF (ANS .EQ. 7 Y 7 .OR. ANS .EQ. 7 y 7 ) GO TO 4002 
IF (ANS .EQ. 7 N 7 .OR. ANS .EQ. 7 n 7 ) GO TO 555 
GO TO 520 



C 

C FORMAT STATEMENTS 

C 

C 



2801 FORMAT 

2802 FORMAT 

2803 FORMAT 
+ 

2804 FORMAT 
+ 

2805 FORMAT 
+ 

2811 FORMAT 

2806 FORMAT 
+ 

2807 FORMAT 
+ 

2808 FORMAT 
+ 

2809 FORMAT 

2810 FORMAT 

2813 FORMAT 

2815 FORMAT 
+ 

2816 FORMAT 

2817 FORMAT 
+ 

2818 FORMAT 
+ 

2819 FORMAT 
+ 

2840 FORMAT 

3000 FORMAT 
+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4 - 

3001 FORMAT 

3003 FORMAT 

3004 FORMAT 

3005 FORMAT 
+ 

+ 



(BN, Al) 

(BN, F15.0) 

(/' Enter the height of the spine in 
, millimeters : 7 , 2X, \) 

(/' Enter the diameter of the spine in 
, millimeters : 7 , 2X, \) 

(/' Enter the thermal conductivity of the 
, spine material ') 

(' in watts per meter per deg C:',2X,\) 

(/ 7 Enter the saturation temperature of the 
',/, 7 fluid in deg C: 7 ,2X,\) 

(/ 7 Enter the temperature of base surface 
, in deg C: 7 , 2X, \) 

(/' The computed spine base temperature is 
, (deg C) : 7 ,5X,\) 

(F5.1) 

(/ 7 The heat dissipated by the spine is 
, (watts) : 7 , 8X, \) 

(/' The revised spine height is (mm) : 7 , 20X, \) 
(/' Percentage of spine in natural convection 
, is : 7 , 4X, \) 

(F9.1) 

(/' Percentage of spine in nucleate boiling 
, is : 7 , 6X, \) 

(/ 7 Percentage of spine in transition boiling 
, is : 7 , 4X, \) 

(/ 7 The spine diameter (user input) is 
, (mm) : 7 ,13X,\) 

(BN, Al) 

(////, 5X, 7 You have entered the following 
, parameters 

5X, 'Spine height (mm) : 7 , 8X, F6 . 12 ,/ , 5X, 

'Spine diameter (mm) : 7 

, 7X, F5 . 1 , / , 5X, 'Thermal conductivity: 7 , 4X, 

F6. 1,/, 5X, 

'Saturation temperature: 7 , 2X, F6 . 1 ,/ , 5X, 

7 Base temperature : 7 , 8X , F6 . 1 , / , 5X , 7 Tip 
temperature: 7 , 

9X, F6 . 1,//// , 5X, 7 Do you want to make any 
, changes? : 7 , 

//,5X, 'Enter "Y" for 7 , 7 yes, "N" for no. 7 ,/) 
(F6.1) 

(BN, Al) 

(BN, II) 

(////, 5X, 'Do you want to enter another data 

t SGt*? , f / f 

5X, 'Enter "Y" for yes, "N" for no. 7 ,/) 



3006 FORMAT 

3007 FORMAT 
+ 

3008 FORMAT 
+ 

+ 



(/,' Enter the desired tip temperature in deg 
, C: 7 , 2X, \) 

(/,' The spine tip temperature is (deg C) : 7 , 
15X , \) 



(////, 5X, 'Do you want an overview of the 
, program? 7 ,/, 

5X, 'Enter "Y" for yes, "N" for no. 7 ,/) 
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555 CALL CLS 
STOP 
END 

C 

c************************end of main program***************** 

c 

c************************************************************ 

c 

c 

C THIS SUBROUTINE DISPLAYS THE CHOICES OF PARAMETERS SO 

C THAT THE USER MAY CORRECT ANY INPUTS 

C 
C 

c 



SUBROUTINE CORRECT 

CALL CLS 
WRITE (IOT, 5000) 

5000 FORMAT (////,5X, ' (1) Spine height' ,/, 5X, ' (2) 

+ , Spine diameter',/, 

+ 5X, , (3) Thermal conductivity' ,/, 5X, ' 

+ , (4) Saturation' 

+ ' temperature' ,/, 5X, ' (5) Desired base 

+ , temperature',/, 

+ 5X,'(6) Spine tip temperature', 

+ /////, 5X, 'Enter the number of the parameter 

+ , you want to ' , 

+ ' change. ' ,/) 



RETURN 

END 

C 

C 

£**************************★********************************* 

c************************************************************ 

c 

c 

C THIS SUBROUTINE INPUTS ANY CORRECTIONS THE USER 

C DECIDES TO MAKE 

C 

c 

c 



SUBROUTINE UPDATE ( INPUT , ELINP , DI AM , CONDTH , TSAT , TBASE , 
+ TTIP) 

CALL CLS 

IF (INPUT .EQ. 1) THEN 
WRITE ( IOT, 5001) 

5001 FORMAT (//, 5X, 'Enter the height of the spine in 

+ millimeters :', 2X, \) 

READ (IN, 5010) ELINP 
ELSEIF (INPUT .EQ. 2) THEN 
WRITE ( IOT , 5002) 

5002 FORMAT (//, 5X, 'Enter the diameter of the spine in' 

+ ' millimeters :', 2X, \) 

READ (IN, 5010) DIAM 
ELSEIF (INPUT .EQ. 3) THEN 
WRITE ( IOT, 5003) 

5003 FORMAT (//, 5X, 'Enter the thermal conductivity of the 
+ spine material in watts per meter per deg C:',2X,\) 

READ (IN, 5010) CONDTH 



c; o 



ELSEIF (INPUT .EQ. 4) THEN 
WRITE (IOT, 5004) 

5004 FORMAT ( // , 5X , ' Enter the saturation temperature of the 

+ , fluid in deg C:',2X,\) 

READ(IN, 5010) TSAT 
ELSEIF (INPUT .EQ. 5) THEN 
WRITE (IOT, 5005) 

5005 FORMAT (//, 5X, 'Enter the temperature of the base 

+ , surface in deg C:',2X,\) 

READ (IN, 5010) TBASE 
ELSEIF (INPUT .EQ. 6) THEN 
WRITE (IOT, 5006) 

5006 FORMAT (//, 5X, 'Enter the spine tip temperature in deg 

+ ,C:',2X, 

+ \) 

READ (IN, 5010) TTIP 
ENDIF 

5010 FORMAT (BN, F15.0) 

RETURN 

END 

C 

C 

C* *********************************************************** 

c* ******************************************** *************** 

c 

c 

C THIS SUBROUTINE DISPLAYS THE PROGRAM OVERVIEW 

C 

c 

c 

SUBROUTINE OPTION 
CALL CLS 



WRITE (IOT, 6000) 

6000 FORMAT (///, 2 OX, '»>»>»> PROGRAM OVERVIEW ««««<' 

+ ,// , 10X, 'This program will predict the 

+ ,performance of a spine', 

+ ' in the' ,/, 8X, 'boiling mode. The boiling 

+ ,liquid for this program model is ',/, 8X, 'water . 

+ , The algorithm divides the spine into 100 

+ ,subfins' ,/,8X, 'and uses a cascade procedure to 

+ , analyze the performance of, /,8X,'the spine. 

+ ,The cascade begins at the fin tip', 

+ ' and concludes at the' ,/, 8X, 'base surface. The 

+ ,user will be asked to enter the',/,8X, 

+ , 'following in SI units: ',//, 10X, 

+ '> Spine height (mm) ' ,/, 10X, '> Spine diameter 

+ , (mm) ' ,/, 10X, '> Thermal conductivity of the 

+ , spine (W/m - deg C) ' , 

+ / , 10X , ' > Saturation temperature (deg C)',/,10X, 

+ '> Desired base temperature (deg C)',/,10X,'> 

+ , Desired tip temperature (deg C)',//,8X, 

+ 'The program may adjust the spine height input 

+ , in order' ,/, 8X, 'to meet the specified base 

+ , temperature. The spine', 

+ /, 8X, 'diameter is not adjusted.') 



WRITE (IOT, 6001) 

6001 FORMAT (//' PRESS ENTER TO CONTINUE ' , 2X,/) 

READ (IN, 6002) RESPON 

6002 FORMAT (BN, Al) 



CALL CLS 
WRITE (IOT, 6003) 

6003 FORMAT (///, 10X, 'The model assumes a saturation 
+ , temperature ' , 

+ ' above 100 deg C.',/, 

+ 8X,'It is suggested that the tip temperature be 

+ , at least 2 deg C' ,/, 8X, 'higher than the 

+ , saturation temperature. A numerical 

+ , instability ',/, 8X, 'may occur', 

+ ' with a low tip temperature excess 10X, 

+ 'If any value is entered incorrectly, the user 

+ ,is given' ,/, 8X, 'the opportunity to change any 

+ ,or even all of the values' ,/,8X, 'before the 

+ , program is run.',/, 

+ 10X, 'The following are the outputs displayed by 

+ ,the program: ',//, 10X, '> Adjusted spine height 

+ , (mm) ',/,10X, 

+ '> Spine diameter (mm) ' ,/, 10X, '> Heat flow 

+ , (Watts)',/, 

+ 10X, '> Predicted base temperature (deg C)',/, 

+ 10X,'> Adjusted tip temperature' ,/, 10X, 

+ '> Percent of spine in natural convection' ,/, 10X, 

+ '> Percent of spine in nucleate boiling ',/,10X, 

+ '> Percent of spine in transition boiling') 

WRITE (IOT, 6001) 

READ (IN, 6002) RESPON 

CALL CLS 
WRITE (IOT, 6004) 

6004 FORMAT (///, 1 OX, 'All values are rounded off to the', 

+ ' tenth decimal place. ',/, 10X, 'At the 

+ , conclusion' , 

+ 'of the program, the user is given the option 

+ , to' ,/, 8X, 'either enter another data set or to 

+ , terminate. ' ) 

WRITE (IOT, 6001) 

READ (IN, 6002) RESPON 

CALL CLS 

RETURN 

END 

c**********************end of subroutines******************** 
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